Random and Correlated Roughening in Slow Fracture by Damage Nucleation 
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We address the role of the nature of material disorder in determining the roughness of cracks 
which grow by damage nucleation and coalescence ahead of the crack tip. We highlight the role of 
quenched and annealed disorders in relation to the length scales d and associated with the disorder 
and the damage nucleation respectively. In two related models, one with quenched disorder in which 
d ~ ^c, the other with annealed disorder in which d <g ^c, we find qualitatively different roughening 
properties for the resulting cracks in 2-dimensions. The first model results in random cracks with 
an asymptotic roughening exponent ( « 0.5. The second model shows correlated roughening with 
( ~ 0.66. The reasons for the qualitative difference are rationalized and explained. 



I. INTRODUCTION 

When cracks develop slowly via the nucleation of dam- 
age ahead of the tip, the crack surfaces left behind ap- 
pear to be rough. A question of great interest is what 
is the scaling exponent that characterizes the roughening 
of such surfaces and how to relate the value of the ex- 
ponent to the physical phenomena that govern the crack 
propagation. For cracks surfaces in 2+1 dimensions the 
anisotropy of the fracture experiment results in a number 
of scaling exponents, making the attainment of a satis- 
factory theory quite difficult T]. On the other hand, for 
cracks in quasi 2-dimensional samples, where the result- 
ing surfaces are 1+1 dimensional graphs (rupture Hues), 
the issues are clear at least in the sense that there exists 
one well-defined scaling exponent. This is conveniently 
defined by measuring y{x) where y is the height of the 
graph above the Euclidean coordinate x that defines the 
crack direction and then defining some measure of the 
height fluctuations, for example 



h{r) = {Max{y{x)}^ 



<a:<a:+r 



Min{y{x)} 



(1) 

For self-aSine graphs the scaling exponent C is defined 
via the scaling relation 



(2) 



It is well known that random graphs are consistent with 
C = 0.5, whereas positively (negatively) correlated graphs 
are characterized by C > 0.5 (C < 0.5). Experiments on 
2-dimensional samples tend to report scaling exponents 
in the range C « 0.65 ± 0.04 0, H 0, indicating the 
existence of positive correlations between successive crack 
segments. 

In recent work a model was proposed in 1-1-1 dimen- 
sions for such slow crack propagation via damage nucle- 
ation and coalescence ahead of the crack tip [3, @|- A 
crucial aspect of this model is the existence of a typical 
length scale ahead of the crack tip where damage nu- 
cleation can take place Q . A physical picture that might 
support such a scenario (though definitely not a unique 
one) is that a small plastic zone of linear dimension 
forms around the crack tip and the relevant damage units 
involved in the process are plastic voids. The idea is that 



since plastic deformation is typically associated with a 
limiting stress level (denoted the "yield stress") a^, the 
purely linear elastic divergent stresses are cut off such 
that they cannot reach a critical level required for the nu- 
cleation of voids at the crack tip. It is claimed that such 
critical levels of stress can be attained approximately 
near the outer boundary of the plastic zone, i.e. the 
elastic-plastic boundary. Whenever a void is nucleated 
in this region, it evolves and eventually coalesces with 
the current crack tip to form a new crack configuration. 
The crack then evolves by successive applications of such 
nucleation and coalescence events. Refs. ^ P| demon- 
strated that the rupture lines generated by this model 
are self-affine rough graphs with a correlated scaling ex- 
ponent C « 0.66. Not only the value of the roughness 
exponent is found to be significantly above the random 
walk exponent ( = 0.5, it also appears close to the mea- 
sured one Q, H, Q . The aim of this paper is to gain a 
deeper understanding of the origin of this result. 

An essential question asked in the context of any such 
fracture growth model is how to represent and incorpo- 
rate the effect of material disorder. This issue is im- 
portant since it asks how small scale features affect large 
scale properties, for example the power-law scaling of Eq. 
^ that indicates a lack of characteristic length scale. To 
our knowledge there had been no systematic study of the 
role of the nature of material disorder in determining the 
roughness of cracks in 1-f 1 dimensions. Our aim here is 
to shed some light on this issue by demonstrating that 
different views of material disorder and the associated 
length scales have a qualitative effect on the scaling prop- 
erties of rupture lines. To this aim we elaborate on the 
type of material disorder adopted in the model described 
briefiy above - referred to below as model A - and present 
a new model - referred to below as model B - that incor- 
porates a different picture of material disorder. In model 
B the disorder is quenched, and the stochasticity associ- 
ated with the material heterogeneities is fixed a-priori in 
space and time. In model A the disorder is "annealed" 
in the sense that the stochasticity depends on the actual 
state of the system. Whenever the disorder has some 
spatial characteristic scale we denote it by d and call it 
the "disorder length" that should be compared to the 
previously introduced length fc- Thus model A is char- 
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acterized by annealed disorder and li <C ^c, while model 
B is characterized by quenched disorder and d ~ ^c- One 
of the points of this paper is that this change in material 
disorder in model B is sufficient to destroy the positive 
correlations between successive crack segments observed 
in model A, changing the universality class of the model 
and ending up with a random graph with an asymptotic 
scaling exponent C « 0.5. 

In Sect. |n] we present in more detail model A and 
recall its results, elaborating on the way in which ma- 
terial disorder is incorporated into the model. We then 
explain the modifications leading to model B. In Sect. IIIII 
we present the new results for model B, compare them 
with model A and clarify the origin of the qualitative dif- 
ferences between them. Sect. II VI offers a summary and 
some concluding remarks. 



II. CRACK PROPAGATION BY DAMAGE 
NUCLEATION AND COALESCENCE 

A. A non-perturbative calculation of the 
linear-elastic stress fields 



Note that even though the initial configuration is that 
of a straight crack, with material disorder the crack 
might deviate from the straight path, attaining an arbi- 
trary rough shape. Solving the bi-Laplace equation with 
boundary conditions on such an arbitrary boundary is 
quite a formidable task. Recently, we have developed a 
general method of solution based on iterated conformal 
maps 0. In this method, one starts with a crack for 
which the conformal map from the exterior of the unit 
circle to the exterior of the crack is known. For exam- 
ple, in our case we start with a long straight crack in 
the form of a mathematical branch-cut, representing the 
common experimental practice of introducing the sample 
with a notch in order to localize the fracture process in 
a controlled way. We can then grow the crack by little 
steps in the desired directions, computing at all times 
the conformal map from the exterior of the unit circle to 
the exterior of the resulting crack. Having the confor- 
mal map makes the exact calculation of the stress field 
straightforward in principle and highly affordable in 
practice. The details of the method and its machine im- 
plementations are described in full detail in 0. 



The mathematical difficulty in developing a theory for 
the morphology of fracture surfaces is the necessity of 
calculating the linear-elastic stress fields for highly non- 
regular crack paths. Typically, the stress conditions near 
the crack tip depend non-linearly on the crack path y{x). 
Formally, one has to solve the bi-Laplace equation for the 
Airy stress potential xi^iU) S 

AAx{x,y)^0, (3) 

in the infinite plane with traction-free boundary condi- 
tions on the crack surfaces 

CTxnis) = (7yn{s) = . (4) 

Here s is the arc-length parametrization of the crack 
shape and the (Ti„(s) denotes the stress acting in the i-th 
direction on a segment whose out-ward normal is the nor- 
mal to the crack face at s. The stress tensor field (Jy is 
derivable form the Airy stress potential x{^^ n) according 
to 

d'^X d'^X d'^X 

The relevant experimental configuration for our purpose 
is that of global mode I fracture in which a system con- 
taining initially a straight crack in the x direction, sub- 
jected to a stress applied in the direction y, perpendicular 
to the crack. At infinity we write the boundary condi- 
tions 

CTxxioo) = ; CTyyioo) = (T°° ; axy{(X>) = , (6) 



B. Model A 

1. Damage nucleation 

We consider a crack evolving under quasi-static condi- 
tions by the nucleation and coalescence of damage ahead 
of the crack tip. These damage elements can be voids 
or micro-cracks. We focus on situations where only one 
damage element nucleates before the process of coales- 
cence. This process is associated with a length scale 
characteristic of the distance of the damage element from 
the tip. A plausible physical picture for such a process 
was proposed in Q . The idea is to identify with 
the size of the plastic zone that develops near the crack 
tip due to the large stresses concentrated there. More 
specifically, it was assumed that the material flows plas- 
tically such as to reduce the stress field near the crack 
tip to a level determined by the yield stress a^. Math- 
ematically, the statement is that the distortional energy 
J2 =^SijSij, with Sij = aij — |Trcr(5,y, satisfies the rela- 
tion [ij 

•h = ol, (7) 

inside the plastic zone. Outside this region, the stress 
field behaves linear-elastically. To find the outer bound- 
ary of the plastic zone, which has a characteristic length 
^c, we use the iterated conformal mapping solution of the 
linear-elastic problem. We calculate the spatial curve for 
which Eq. Q is satisfied when approaching the crack 
tip. This curve defines the elastic-plastic boundary. It 
was further shown |^ 0| that the hydrostatic tension P, 
defined as 



where a°° is assumed to be constant. 



(8) 
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attains a larger value near the elastic-plastic boundary 
than inside the plastic zone. Under the physically plau- 
sible assumption that damage will nucleate in regions 
where P exceeds some threshold value Pc, we expect 
damage to nucleate near this boundary. As was explained 
before, after damage nucleates it evolves such that it co- 
alesces with the tip, generating a new plastic zone under 
the influence of the liner-elastic fields and so on. Note 
that in this physical interpretation the damage elements 
are plastic voids and the coalescence process is assumed 
plastic as well (for example, necking of the ligament be- 
tween the crack tip and the void). In fact, as we are not 
resolving the processes by which the crack tip coalesces 
with the void ahead of it, using the nucleation site only 
as a pointer for the advance of the crack. Therefore, we 
are only interested in the roughness of the crack on scales 
larger than ^c- 



2. Growth rule and Disorder 

Naturally, the precise location of the nucleating dam- 
age may be stochastic due to material disorder. To 
quantify this, assume that nucleation occurs only at lo- 
cations r in which the hydrostatic tension P exceeds 
some threshold value Pc- Given the distribution P{r) 
we consider a probability density function f{P{r) — Pc). 
One has in mind an activation process for the nucleation 
of damage, and this activation is more efficient when 
P{r) — Pc is large. The probability for activation vanishes 
for P — Pc < 0. This activation may be due to stress cor- 
rosion in one case, or due to other mechanism in another 
case, but the important thing to note is that P{r) is a 
long-ranged functional of the history of crack evolution, 
potentially leading to the long range correlations implied 
by ^ > 0.5. The fact that damage nucleates depending 
on the stress field aij through P without reference to any 
pre-determined distribution of disorder implies that the 
model is characterized by annealed disorder. Note that 
the same formulation describes equally well a situation in 
which damage nucleates at points where the material is 
weak, if the random weak points are dense enough such 
that the typical scale d separating them is much smaller 
than ^c and that the distribution of nucleation thresh- 
olds is immaterial, characterized only by Pc. The relation 
d <ti ^c allows us to take the continuum limit to define 
a probability distribution function. Model A was stud- 
ied in Refs. 0, 0| . In the absence of precise knowledge 
of the activation process we adopted reasonable probabil- 
ity distribution functions f{P{r)~Pc) and demonstrated 
that the cracks generated by the model were self-affine 
with C = 0.66 ± 0.03 irrespective of the specific form of 

fim - Pc). 

This model should be contrasted with the more com- 
mon mathematical representation of stochastic growth 
models via a Langevin type equation. In this case a de- 
terministic equation is supplemented with an additive 
noise term whose statistics are independent of the de- 



terministic part. In model A the randomness cannot be 
represented by an additive independent noise. We now 
turn to model B to test the influence of the type of ran- 
domness employed on the roughness of cracks. 



C. Model B 

In model B the crack is still assumed to propagate by 
the nucleation and coalescence of damage ahead of its 
tip. The linear-elastic stress fields are still calculated us- 
ing the powerful method of iterated conformal mapping 
described in Sec. Ill Al The main difference between the 
two models stems from a different way of incorporating 
material disorder into the crack evolution process. In 
model B the disorder is assumed to be quenched, rep- 
resented by an a-priori random distribution of identical 
weak points. The random weak points have a prescribed 
density such that the average distance between them is 
d. Physically, the weak points can be realized by den- 
sity fluctuations in an otherwise homogeneous material 
or by small particles that have a lower breaking thresh- 
old than the matrix in which they are embedded, but 
do not change signiflcantly the elastic properties of the 
system. As in model A, the damage nucleation process 
near the crack tip is characterized by a length scale ^c- 
A second point of departure from model A is that in this 
model we assume d ~ This relation can be realized 
in different physical situations. For example, can be 
still identified (as in model A) with the linear dimension 
of the plastic zone, where the independent scale d just 
happens to be of the order of magnitude; in that case the 
damage elements can still be plastic voids. Alternatively, 
if plastic process are not dominant, one could imagine 
the crack pinned to a weak point until a micro-crack nu- 
cleates at another weak point to propagate the crack by 
coalescence. In this interpretation is by definition of 
the order of d. 

To complete the model we need a growth rule. As was 
mentioned before, the weak points are assumed identical 
in the sense that they have the same breaking threshold 
that is sig niflcantly smaller than the ordinary material 
points . Since the disorder in this model is quenched, 
we can define a deterministic growth rule stating that 
the crack advances to the weak point where the hydro- 
static tension P is maximal. Note that even though the 
weak points are spread randomly and independently of 
P, the selection of a weak point to be a pointer for the 
next crack growth depends crucially on its spatial prox- 
imity to the maximal hydrostatic tension P; the closest 
weak point to the maximal hydrostatic tension P is most 
likely to be chosen at each growth step. It is worthwhile 
mentioning that if one could fix and decrease d such 
that d/S^c (a limit that is not realized in our model 
where d ~ ^c), one would obtain the deterministic limit 
of the model since the maximal hydrostatic tension P 
would almost inevitably coincide with a weak point and 
the crack would advance almost always to the point of 
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maximal hydrostatic tension. Thus, the ratio d/^c is a 
measure of the width of the statistical distribution in this 
model. In the next section we analyze the new model and 
compare its results to the results of mode A. 



III. RESULTS AND DISCUSSION 

We have simulated model B and obtained several crack 
realizations each of about 500 growth steps. An example 
of a resulting crack is shown in Fig. ^ 
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FIG. 1: A typical crack of 500 growth steps in a.u. Note the 
difference in scales between the ordinate and the abscissa. 
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FIG. 2: A log-log plot of the height fluctuations h{r) of Eq. 
0as a function of the scale r. The straight line corresponds 
to C = 0.57. 



Common (and non-trivial) properties of models 
A and B 



The first thing to point out is that both models A 
and B deviate from many models common in the litera- 
ture in a way that might affect significantly the scaling 
properties. In both models the crack propagates via the 
nucleation and coalescence of damage at a finite distance 

ahead of its tip, essentially under the influence of a 
linear-elastic stress field. In ideal linear elasticity, the 
stress tensor field ct^ attains the following asymptotic 
form approaching the crack tip 



We have measured the roughness of the cracks in 
the model with both the variable bandwidth max-min 
method of Eq. ^ and the variable bandwidth RMS 
method Ts'l . In order to avoid strong finite size effects we 
have used the results of Ref. d3] to calibrate our results 
for the different measurement methods. This procedure 
turned up to be consistent in the sense that the variance 
in the results obtained by the different methods pointed 
to a single well-defined exponent according the finite size 
effects predicted in . Finally, we averaged over differ- 
ent realizations, obtaining ( = 0.53 ± 0.03. An example 
a single roughness measurement is shown in Fig. [3 

This result indicates that the cracks in model B ex- 
hibit random roughness as the roughness exponent is not 
significantly different from the random walk exponent 
C = 0.5. This result is qualitatively different from the 
results of model A in which ( — 0.66±0.03. At this point 
we must conclude that even though the two models share 
many features and on the face apparently one should not 
expect dramatically different scaling properties, the two 
models belong to different universality classes. Model A 
exhibits correlated roughening, while model B exhibits 
random roughening. We turn now to further clarifying 
the origin of the qualitatively different results. 



cTy (r, 6) = 



(9) 



where (r, 9) is a polar coordinates system located at the 
crack tip, and are known universal functions and 
Kj and Kn are the stress intensity factors correspond- 
ing to opening (mode I) and shearing (mode II) stresses 
Eq. ||2J) describes well the stress fields on a scale r 
relative to the crack tip that is much smaller than any 
other length scale in the problem. Naturally, we cannot 
expect this formula to be precise for r of the order of ^c- 
In particular, the hydrostatic tension P{r) should be sen- 
sitive to significant corrections to the ideal law Eq. 1^. 
We can expect that the stress fields on a scale away 
from the tip in both models is described by Eq. ^ with 
Kj ^ Kji plus additional terms. To verify this expecta- 
tion we calculated the hydrostatic tension P on an arc a 
distance from the tip and fitted to the form 



Pi9) 



Ki 



Kii . 

sm , 

Ki \2 



(10) 



predicted by Eq. (|UJ) The results support our ex- 

pectation, showing that Kj ^ Kjj and an additional 
contribution of 10-15% from other non- universal terms 
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FIG. 3: The hydrostatic tension P (see Eq. as a func- 
tion of at a fixed radius r = (solid line) and a fit with 
Ki = 125 {a.u.), ^ = -0.51 (dashed line). The present 
data was taken from a crack after 300 growth steps, but it is 
representative. 



at a distance away from the tip for both models. An 
example is shown in Fig. |31 

The conclusion is that both models are in the non- 
linear regime (in terms of Kji{y{x)}). In addition one 
cannot neglect the influence of contributions on top of 
Eq. We believe that this is an important aspect in 
the success of model A in reproducing a correlated rough- 
ening, quantitatively close to the experimental observa- 
tions. It also explains why models which used the linear 
approximation of a straight crack plus a perturbed Kn 
achieved scaling exponents different from those observed 
in experiments 15,16]. 

On the other hand, this cannot be the only factor ex- 
plaining the results of model A since this feature is com- 
mon to model B which fails to produce correlated rough- 
ening. 



B. The difference between the models stems from 
the difference in randomness 

We propose that the differences in randomness are re- 
sponsible for the different universality classes of model 
A and B. Recall that the quenched disorder length d in 
model B is of the order of the nucleation length ^c- As a 
result, there are only few weak points available for dam- 
age nucleation on a scale ^c- In this situation the crack 
"selects" a damage nucleation site out of a small number 
of possibilities and in fact there is a sizeable probability of 
having a growth step chosen uncorrelated with the deter- 
ministic field P that carries the history of the evolution. 
We thus expect this randomization effect of uncorrelated 
growth steps to accumulate gradually and reduce the cor- 
related roughening exponent observed in model A. The 



situation is fundamentally different in model A where the 
probability of having completely random growth steps is 
small. To support this explanation we have measured 
the roughness exponent C (see Eq. |2l as a function of the 
crack length in terms of growth steps. In Fig. 0]we show 
the dependence of C on the number of growth steps. 
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FIG. 4: The average roughness exponent as a function of 
the number of growth steps N. 



It is observed that indeed roughness exponent ^ de- 
creases monotonically from C ~ 0.625 to C 0.53, pre- 
sumably approaching asymptotically C « 0.5. This find- 
ing is consistent with the suggested explanation. 



IV. SUMMARY 

We have studied the role of disorder in generating cor- 
related roughening exponent in cracks of 1 -I- 1 dimen- 
sions. The main conclusion is that adding additive ma- 
terial disorder to linear elasticity is not sufficient to gen- 
erate correlated crack graphs with exponents larger than 
0.5. This is due to the destructive events or uncorrelated 
steps which accumulate and gradually produce a random 
graph. The generation of correlated graphs (C > 0.5) 
is due to the correlations between the deterministic field 
(the hydrostatic pressure) and the pdf carrying the ran- 
domness, like the annealed disorder of model A. We re- 
iterate the result concerning the magnitude of the stress 
intensity factors, and specifically Kjj. The measurement 
of Kjj and higher order terms of the pressure field tells 
us that linear approximations of Kjj, or the principle 
of local symmetry, do not represent properly the stress 
field at the vicinity of away from the rough crack tip. 
Therefore, one can view this feature as an additional test 
for models of crack growth which produce the observed 
correlated roughness. 
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